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Abstract 

A new approach to high energy evolution based on a linear equation for QCD generating 
functional is developed. This approach opens a possibility for systematic study of correlations 
inside targets, and, in particular, inside realistic nuclei. 

Our results are presented as three new equations. The first one is a linear equation for 
QCD generating functional (and for scattering amplitude) that sums the 'fan' diagrams. For 
the amplitude this equation is equivalent to the non-linear Balitsky-Kovchegov equation. 

The second equation is a generalization of the Balitsky-Kovchegov non-linear equation to 
interactions with realistic nuclei. It includes a new correlation parameter which incorporates, 
in a model dependent way, correlations inside the nuclei. 

The third equation is a non - linear equation for QCD generating functional (and for scat- 
tering amplitude) that in addition to the 'fan' diagrams sums the Glauber-Mueller multiple 
rescatterings. 
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1 Introduction 



High density QCD [H El El IH E] which is a theory of Color Glass Condensate deals with parton 
systems with large gluon occupation numbers. It has entered a new phase of its development: a 
direct comparison with the experimental data. A considerable success (HI 13 El EO CHI EH EE21 UH1 
HH HH] has been reached in description of new precise data on deep inelastic scattering [16^ as well 
as in understanding of general features of hadron production in ion-ion collision J7| . 

Most of the applications of hdQCD are based on or related to a nonlinear evolution equation 
derived for high density QCD [TJ El CHI UH] • The equation sums 'fan' diagrams, which is a subset 
of the semi-enhanced diagrams. In order to obtain a closed form equation it is usually assumed 
that color dipoles produced by the evolution interact independently [TJ [THl - It was noticed in Ref. 
[TT] that the nonlinear evolution equation based on the above assumption faces problems when 
applied to realistic nuclei. It was realized that in addition to multiple rescatterings on nucleons 
inside a nucleus we have to include in the analysis multiple rescatterings inside a single nucleon. 
This means a need to account for nucleus correlations. It is our motivation for this research to 
develop a systematic approach which would allow us to introduce correlations inside a target. 

We have found that a convenient language to address the problem of target correlations is 
using the QCD generating functional, for which a linear evolution equation can be written. The 
generating functional was introduced by Mueller in Ref. [2U] and then used by Kovchegov in his 
derivation of the non-linear evolution equation for interaction amplitude ^H]- in the present paper, 
we show that this functional (and the amplitude) obeys a linear evolution equation involving 
functional derivatives with respect to initial conditions. For the linear equation we propose a 
systematic method to include target correlations. In general, the linear equation for the amplitude 
with target correlations cannot be transformed to a non-linear one. For a specific choice of 
correlations we are able to reduce the linear equation to a non-linear form similar to the Balitsky- 
Kovchegov (BK) equation. In fact, our new equation is a model dependent generalization of the 
BK equation for targets with correlations. It is important to stress that we consider the very same 
functional as of Mueller. The evolution of the dipole wave function is still based on independent 
production of dipoles. Only correlations for the dipole-target interaction are introduced. 

Let us discuss now realistic nuclei and reasons why we believe target correlations are essential 
there. In principal, the BK equation is correct for hadron targets as well as for heavy nuclei. 
For the latter the large A limit is usually assumed which adds an additional justification to the 
equation validity. In this case a nucleus is viewed as a very dense system of nucleons. We are 
going to relax this assumption. In fact the realistic nuclei are not very dense but rather dilute 
|1U1 The diluteness parameter that governs the interaction with nuclei could be expressed in 
the form: 



where ttR 2 n is the effective area in which gluons of a nucleon are distributed. 5^(6) stands for a 
nucleus profile function, which can be taken as the Wood-Saxon form factor of the nucleus. 

The parameter ka oc and , in principle, is large for heavy nuclei. It turns out, however, 
that the gluons are distributed inside a nucleon within a rather small area with radius of the order 
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of 0.3 — 0.4/m ^01 E3, which is much smaller compared to the proton radius. For such small 
area the value of ka even for the heaviest nuclei is smaller than unity. Therefore, in spite of the 
pure theoretical interest to the very heavy nuclei we have to be very careful applying the large A 
limit to the realistic nuclei. It should be stressed that even for more 'standard' approach with the 
radius of gluon distribution of the order of 0.6 fm the value of ka reaches at best 3 for heaviest 
nuclei. In this case the 1/ka corrections are of the order 35 % and should be accounted for when 
confronting with experimental data. We remind that the most interesting effects in ion-ion and 
hadron-ion collisions such as Cronin enhancement and A/p art -scaling J7j are of the level about 
25%. 

In this paper we consider realistic nuclei without the large A assumption. We argue that in 
this case a nonlinear evolution is still a valid tool. However the equation gets modified compared 
to the one of Ref. jTHI with which our discussion will have a strong overlap. For realistic (dilute) 
nuclei the assumption about independent dipole interactions should be relaxed. Within some 
model we will be able to include correlations inside nucleus such that only one correlation factor 
is introduced. This factor is independent of the number of dipoles interacting with the nucleus. 

The problem which we are going to deal with is very simply illustrated on the following toy 
example 1 . Suppose we have a probe with two dipoles only which are close to each other on the 
inter-nucleon distance scale T4 ~ A 1 / 3 . Suppose the effective interaction radius of a nucleon is 
i?7v- Then the probability for one dipole to interact with the nucleus is AR 2 N jT\. What will 
be the probability for two dipoles to interact with the nucleus? For the case of independent 
scattering it is (AR 2 N /T^) 2 . It is obvious, however, that this is not what actually happens. Since 
the dipoles are close to each other they either both hit the nucleon or both miss it. The real 
probability will be rather AR 2 N jT\. It follows from this discussion that the equation based on 
independent interactions underestimates shadowing corrections in nuclei. In the above example 
K A = AR 2 N /T 2 . 

If schematically the equation for the amplitude of Ref. ^H] has the form 



dN K 



N — N 2 



1.2) 



then the equation which we propose for the dipole-nucleus amplitude can be written as 



dN A „ 

= Ker 

dY 



K A 



;i.3) 



Eq. (jl.3|) can be certainly brought to the form of Eq. (jl.2|) provided a corresponding rescaling of 
the initial conditions is done. When ka is large Eq. ()1.3|) reduces to Eq. (jl.2|) . In the opposite 
limit ka is small, Na = i^aNn with being a solution of Eq. ()1.2|) for a nucleon target. 

In the final part of this paper we generalize the approach based on QCD generating functional 
to include Glauber-Mueller multiple rescatterings. In order to achieve our goal we construct a QCD 
effective vertex for one dipole splitting to many. In Pomeron language, this vertex generalizes the 
famous triple Pomeron vertex to the case of multiple splitting. We define a new QCD generating 

*We are thankful to Alex Kovner who understood our problem and brought in this example. 



3 



functional which obeys a new linear evolution equation. Without target correlations this equation 
can be reformulated as a new non-linear equation. In addition to the 'fan' diagrams, the functional 
sums all types of semi-enhanced diagrams. We argue that in the saturation domain all semi- 
enhanced diagrams are of equal importance. This is why they have to be summed up in order to 
get a reliable description of the saturation domain. We find that the inclusion of the Glauber- 
Mueller multiple rescatterings leads to a significant modification of results based on 'fan' diagrams 
(BK equation) only. 

The structure of the paper is following. In the next section (2) we derive a linear evolution 
equation for generating functional and the interaction amplitude including target correlations. 
Section 3 is devoted to realistic nuclei. A model for nucleus correlations is proposed and a gener- 
alization of the BK equation is given. In section 4 we present a generalization of the functional 
approach to Glauber-Mueller multiple rescatterings. The concluding section 5 contains discussion 
of the results obtained. 

2 Linear Equation 

Though we usually consider a nonlinear evolution for the scattering amplitude it is possible to 
formulate the very same problem in terms of a linear functional equation. A linear formulation 
appears to be more suitable for treatment of correlations inside a target. We first start from a toy 
model which simplify the whole problem to such an extent that it could be solved analytically. 
Then we proceed to the complete analysis of the QCD generating functional. 

2.1 A toy model 

In pQCD the Pomeron is described by the linear BFKL equation [21] which is integro-differential 
equation. Instead we will consider a phenomenological Pomeron given by 



This equation reproduces the main property of the real BFKL equation, namely, the power-like 
increase of the amplitude as a function of energy. So far we neglect the fact that the actual degrees 
of freedom are color dipoles having different sizes. We will come back to the real QCD in the next 
subsection. 

Let P n {y) be a probability to find n-dipoles with rapidity y in the wave function of the fastest 
(parent) dipole moving with rapidity Y > y. We introduce in a very simple way a probability for 
one dipole to decay into two: 



In f!2.5|) we have assumed that the triple Pomeron vertex is just u ~ & s . For P n (y) we can easily 
write down a reccurent equation (see Fig. Q ) 




(2.4) 




(2.5) 




y 



(2.6) 
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Figure 1: The equation for the probability to find n-dipoles in one fast dipole. The wave line 
shows a single dipole. 

The first term in the r.h.s. can be viewed as a probability of the dipole annihilation in the rapidity 
range (y y — dy). The second is a probability to create one extra dipole. Note the sign minus 
in front of dP n (y)/dy. It appears due to our choice of the rapidity evolution which starts at the 
largest rapidity y = Y of the fastest dipole and then decreases. 

It is useful to introduce the generating function |221 1201 122] 

Z(y,u) = Pn(y)u n , (2.7) 

n 

At the rapidity y = Y there is only one fastest dipole, that is P±(y = Y) = 1 while P n> i(y = 
Y) — 0. This is the initial condition for the generating function: 

Z(y = Y) = u. 

At u = 1 

Z{y, u = 1) = 1 (2.8) 

which follows from the physical meaning of P n as probability. 

Eq. ()2.6)1 can be rewritten as the equation in partial derivatives for the generating function 
Z(y,u): 

dZ(y, u) dZ(y, u) 

- q = -u u{l - u) (2.9) 

ay ou 

The function Z obeying Eq. ()2.9|) with the above given initial condition 



u e -v {Y~y) 

Z(Y-y,u) = ; ^ — . r. (2.10) 

It is important to observe that Eq. ()2.9|) can be rewritten in the nonlinear form. The general 
form of this equation is such that the solution to this equation is Z{u(y))\ when Z{u(y)) is 
substituted into Eq. (J2.9j) . the derivatives dZj du cancel and we are left with an ordinary differential 
equation for u(y). Using initial condition at y = Y we can write 

- jj\v=Y = -wo (Z - Z 2 )\ y=Y . (2.11) 
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Due to implicit dependence on y this form will be preserved at any value of y. Therefore, the 
linear equation (see Eq. (J2.9)0 can be rewritten in the non-linear form: 

dZ(y, u) / 9 \ , 

= -u Q (Z - Z 2 ) (2.12) 



dy 

Eq. (j2.12|) coincides exactly (within the same approximation neglecting the dipole sizes) with the 
equation for the generating functional for the light cone dipole wave function obtained by Mueller 
20\. Eq. ()2.9j) can be viewed as a linearized version of the non-linear Eq. (|2.11j) . 

To find the interaction amplitude we can follow the procedure suggested in Ref.|19j. namely, 

00 1 <Y l 7 

N A (y,b t ) = ImA e \y,b t ) = - £ ^T^U^t). (2.13) 

^ n! [du) n 

Here 7(6*) is minus imaginary part of the elastic amplitude for a single dipole interaction with a 
nucleus at fixed impact parameter b t . The most important assumption made in Eq. ()2.13j) is in 
the independent interaction of n dipoles which is expressed via 7" factor. In the present work we 
will relax this assumption introducing target correlations. Note that the amplitude can be 
found from the following relation 

N A {y, b t ) = 1 - Z(y, j(bt) + 1) . (2.14) 

Let us now come back to Eq. ()2.9j) from which we are going to derive a corresponding equation 
for the amplitude N A . In order to get a closed equation, the linear Eq. (|2.9|) is differentiated n 
times with respect to u and then set u — 1: 

- i-^U = n ZW U + Z<»-')|„ rf (2.15) 

u dy 2 

where = d n Z/du". 

The equation for the amplitude can be easily obtained from Eq. (|2.15|) . The following four 
operations are performed, i) Since Z is a function of Y — y, we make a substitution dZ n /dy = 
— dZ n /dY. Hence the amplitude evolution will be with respect to initial rapidity of the fastest 
dipole Y. (ii) At the rapidity of the target y — yo initial conditions for the interaction of the 'wee' 
dipoles with the target should be specified, (iii) For independent interactions with the target (ii) 
means we multiply Eq. ()2.15|) by the factor 7™. (iv) Divide both sides by n\ and then sum over n. 
The result is 

1 dN A (Y) dN A (Y, 7) 2 dN A (Y, 7) , , 

uj dY 7 <9 7 7 <9 7 1 ' ' 

Eq. (J2.16|) can be rewritten in the non-linear form using the initial condition N A {Y = y ) = — 7: 

- = N A (Y) - N A (Y) (2.17) 

OJq oY 

Eq. 1)2.17)1 is just the BK non-linear equation ^Hl El] in our simplified model. 
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Let us now discuss a way the target correlations can be introduced in the analysis. For the 
linear evolution they appear when we specify in the initial conditions how n 'wee' dipoles interact 
with the target. An appropriate modification is to introduce correlation factors C n such that 

00 1 r) n 7 

N A (y,b) = -E^U7"C, (2.18) 

The coefficients C n are to be found from a nucleus model. The linear equation for the amplitude 
now reads 

ujq oY a 7 y o 7 J a 7 

Where function F is equal to F(n) = C n /C n _i and describes correlations in the interactions of 
n dipoles. In Eq. (|2.19J) . F(j j-) should be understood as operator series expansion. In case F is 
constant the linear equation (|2.19j) reads 

i!W = ,«i tf ,.!f! (2.20) 
UJq oY #7 c/7 

At Y — y the initial condition is N(Y = y ) = — C 7. It is used to obtain from Eq. ()2.20|) the 
non-linear equation: 

- W = Na{Y) ~ K * N a(Y)- (2-21) 

UJq oY 

with Ka = F/Cq. In fully non-correlated case F — Cq — 1, Eq. (j2.21|) reduces to the BK 
equation for the above toy model. In the next section we will consider realistic nuclei for which 
correlations can be modeled. We will show that within certain assumptions the correlation factor 
Ka is indeed a constant but different from unity. It is important to stress that Eq. (J2.2(Jj) can 
be brought to a nonlinear form for constant F only. Otherwise Na is not a function of a single 
variable j(Y) and Eq. (J2.2(J|) is a linear equation involving high orders of partial derivatives. 

The success of the simple toy model considered above leaves us with two hopes for real QCD 
dynamics. First, we believe in a possibility to find a linear functional equation for the generating 
functional. A linear equation is more transparent and more feasible to solve compared to the 
non-linear one. Second, we believe that Eq. (j2.21|) (generalized to include BFKL kernel) is the 
correct generalization of the BK non-linear equation applicable to targets with correlations such 
as realistic nuclei. 

It turns out that a generalization of a linear functional approach to evolution kernels but 
simplified to a double log accuracy is rather straightforward [2H]- For the complete BFKL kernel, 
however, this is a nontrivial problem with which we will deal below. 

2.2 QCD generating functional for 'fan' diagrams 

In this subsection we consider a linear generating functional for 'fan' diagrams in QCD. We include 
both the LO BFKL evolution of the dipole wave function and QCD triple Pomeron vertex (in large 
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N c approximation). Then we essentially repeat the same derivation as for the toy model above. 
For simplicity of presentation we will not treat the impact parameter (b t ) dependencies in a correct 
way which would imply taking care of all relevant shifts in b t . Instead, we will consider b t just as 
an external parameter which often will not be written explicitly. It should be of no problem to 
reconstruct a correct bt dependencies in the final results. 

As in Ref. j2D| we introduce a generating functional Z 

/n 
P n (Y - y, r, b t - n , r 2 ,...,n...,r n ) J[ u{n) d 2 n (2.22) 

71=1 1=1 

where u(fi) = u% is an arbitrary function of r^. P n stands for a probability density to find n dipoles 
with transverse sizes r\, r 2 , . . . , t% ... r n and rapidity y in the wave function of the fast moving 
dipole of the size r and rapidity Y > y. Note that contrary to Refs. [2H EI| we define P n as a 
dimensionful quantity. The functional ()2.22|) obeys two conditions: 

• At y = Y, P n =i = S 2 (f — r*i) while P n >i = 0. In other words, 

Z{Y - y = 0, r, b t ; [u]) = u{r) . (2.23) 



• At u = 1 

Z(Y -y,r, b t ;[u=l}) = 1. (2.24) 

This equation follows from the physical meaning of the functional: sum over all probabilities 
equals one. 

Contrary to the toy model, in QCD the probability to survive for a dipole of size r; is not constant 
but equal to (see Ref. [20] for details) 



ol t 2 

a s u(ri) = — i / — ' d 2 r' = a s \n(r 2 /p 2 ) 

2tt J P [Ti — r) z r u 



with p being some infrared cutoff and as = asN c /ir. The probability for a dipole of the size 
f\ + r 2 to decay into two with the sizes T\ and r 2 is equal to 



a s (ft + r 2 ) 2 
2 ix r\r\ 



Using these two probabilities we can write down the equation for P„: 
dP n (Y - y, r, b t ; n, r 2 . . . n . . . r n ) 



a s dy 



^( r i) p n (Y -y,r,bt 1 r x ,r 2 ...r i ... r n ) 



i=i 



{n + r n ) 



2 



+ E ro ^ 22 P «-! ( F - r ' h ^ ri ' r2 • • • + r n )... tv-i) (2.25) 



S 



Let us now introduce two operator vertices 

V^r, [u]) = a s uo{r) u{r) — ^ (2.26) 

and 

fy v 2 () 

V^ 2 (r, r', [«]) = -i ^ u(0 n(f - ?) -— - . (2.27) 

2 7r r 2 (r — r) 1 ou(r) 

The functional derivatives with respect to it(r) play a role of an annihilation operator for a dipole 
of the size r. The multiplication by u(r) corresponds to a creation operator for this dipole. 

Multiplying Eq. ()2.25|) by the product Yli =1 Ui and integrating over all we obtain the following 
linear equation for the generating functional: 



dV V^i(r\ [u])Z + / dV dV Vi_> 2 (r", r', [u]) Z . (2.28) 



<9Z 

Like in the case of previously considered toy model, Z is a function of a single variable u(y): if 
Z(u(y)) is substituted into Eq. ()2.28|) we find in the l.h.s. dZ/dy = J d 2 r (5Z/Su(r)) (du(r)/dy). 
The integrals / d 2 r 5Z/5u(r) cancel on both sides of the equation and we are left with an ordinary 
differential equation for u(y). With the help of the initial conditions (Eq. (|2.23|0 . Eq. (|2.28j) can 
be rewritten in the non-linear form reproducing the same equation for Z as in Ref. |2*U] : 

d Z 

-u(r) Z{Y - y,r,b t ;[u}) (2.29) 



a s dy 



+ J ^ ~n jfl ^ Z(Y-y, r>, b t ; [„]) Z (Y — y, (f - r*), b t ; [«]) . 
The amplitude A 7 ^ is defined similarly to Eq. (|2.13|) (see Ref. |H]): 



oo -r „ n / 

iV A (y,r, & t ) = • ■ E - / II ^ 2 ^(n;6t) 

n=l n - J i=l \ 



5ui 



Z{Y -y,r,b t ;[u})\ u=1 . (2.30) 



As above C n describe correlations in the n-dipole target interaction. In order to derive an ex- 
pression for the amplitude iV^ we apply to both sides of Eq. ()2.29|) n functional derivatives with 
respect to u and then set it — 1. Using the notation Z^ = ri?=i t~ Z (Y — y, r, b t ; [it]) 



wc 



Sui 

obtain 

(2 ' 31) 

The evolution equation for the amplitude is obtained through the same steps as in the toy model 
above. We divide Eq. ()2.3H) by n\ and multiply by the product — C n Iir=i 7( r i; h)- Then we 
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integrate over all r, and sum over n. Finally we consider the evolution with respect to Y instead 
of y which implies a change of sign in front of the derivative: 



dN A {Y- [7]) 
a s d Y 



-J d 2 r' V^ir', [ 7 (r')]) N A (Y; [7]) (2.32) 



+ F ( [ d 2 f~f(f) -^-] f d 2 r" d 2 r' V^ 2 (r', r", [ 7 (r')]) A^F; [7]) 



5 7 (r) 

Eq. ()2.32|) is linear but in functional derivatives and it is not easy to develop direct methods to 
solve it. In case of function F being a constant, Eq. (|2.32j) can be easily reduced to a non-linear 
equation. Using the initial condition N A (Y = y ) = —Co 7 we obtain the following nonlinear 
equation: 

dN A {r, y; b t ) f d 2 r' r 2 

^7 = a s x — 7 r— — 2 x (2.33) 

oY J P 2 n r 2 (r — r) 2 

[2N A (r', y; b t ) - N A (r, y; b t ) - K A N A (r', y; b t ) N A (f - 7, y; b t )} , 

with K A = F/Cq. If we consider dipole interaction as fully uncorrelated F = Cq = 1 then 
Eq. ()2.33|) reduces to the BK equation. We will find below that within some assumptions, cor- 
relations in realistic nuclei can be described by a constant F, but with F 7^ Cq. In this case 
Eq. ()2.33|) is a generalization of the BK equation to realistic nuclei with correlations. 



3 Dipole-nucleus interactions with correlations 
3.1 Nucleus correlations 

Let us discuss the target correlations starting with dipole-nucleus interactions at sufficiently low 
energy (y = ln(l/x )). It should be stressed that the energy is supposed to be not too low such 
that high energy QCD is applicable. On the other hand, we assume that energy is low enough 
such that a single dipole interacts with one nucleon only. We neglect all Glauber rescatterings at 
this energy. 

Consider first the Fig. El which presents in details how a single dipole interacts with a nucleus. 
Schematically the same is displayed in (Fig. Ola). Its contribution is proportional to K A (b t ): 

J d 2 b'N N (r, y ; b t - b' t ) S A (b' t ) = S A (b t ) J d\" N N (r, y o] 6 t ") « N N (r, y ; b t = 0) K A (b t ) . 

(3.34) 

Here stands for the imaginary part of the dipole-nucleon elastic amplitude. can, in 

principal, include multiple dipole-nucleon rescatterings. As it is clear from Eq. ()3.34|) and Fig. El 
that bt is a distance from the nucleus center to a given nucleon. 
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In Eq. (|3.34|) the first equality is due to the assumption, standard for the Glauber theory, that 
the size of the nucleon as well as the radius of the dipole-nucleon interaction is much smaller than 
the size of the nucleus. The second equality is model dependent. We have assumed a cylindrical 
nucleon that is Nw(r, y ; b t ) = N^(r, y Q ; b t = 0) Q(Rn ~ h). The parameter ka(^) (Eq- Q1-1J) ) 
measures the strength of the dipole nucleus interaction. 

Consider now diagrams of Fig. E3-b. For very large nuclei the diagram of Fig. Ol-b-2 dominates 
since it is proportional to n A (b t ) while the diagram of Fig. OJ-b-1 is of the order K A (b t ) only. In 
the case when HA^h) is not large the diagram of Fig. 01-b-l should be also taken into account. Its 
contribution has the form: 

- J d 2 b' t (N N (r,y ;b t -b' t )y S A (b' t ) ~ (N N (r, y ; b t = 0)) 2 K A {b t ) . (3.35) 

Fig. EJb should be viewed as interaction of two dipoles with the nucleus. Fig. EJJ-b-l is a correlated 
interaction with one and the same nucleon while Fig. 01 b- 2 is independent interaction with different 
nucleons. In Eq. (J3.35J) we have neglected the difference in the dipole sizes, which is not important 
for this discussion. 

If we proceed further for higher order terms we can find a general expression for n-th term 
in the expansion. For three dipoles we have a term proportional to ka when all three interact 
with the same nucleon, a term proportional to 4 when two dipoles are correlated, and a term 
proportional to 4 when all interact independently. The sum of these terms will be 

N 3 N (k a + 24 + 4). 

The factor 2 in front of 4 ls due to combinatorial counting. The n-th term in the expansion has 
the form: 

(_!)» KA (i + KA )«-i jv«. (3.36) 

There is no time ordering associated with the factor 1/n! which usually appears in the Galuber 
expansion. (— 1)" stems from the fact that the dipole-nucleon amplitude at y = ln(l/a;o) ^> 1 
is pure imaginary. 

From Eq. (|3.36|) we can read off the correlation coefficients C n introduced above: C n = k,a (1 + 
^a)" 1 - Hence the correlation function F = (1 + k,a) and the correlation parameter Ka = 
(1 + k>a)/ka- The correlations appeared when the assumption about independent interactions 
was relaxed. In spite of the introduction of the correlation parameter, Eq. ()3.36|) still allow 
interpretation in terms of independent dipole-nucleus interactions. The strength of this interaction 
is (1 + ka) Nn- In fact the dipoles interact 'almost' independently and only the overall factor 
+ &a) incorporates the correlations inside the nucleus. The most important fact is that 
we were able to introduce only one correlation parameter independent of the number of dipoles 
participating in the interaction. This success can be traced back to the model assumption about 
the cylindricity of nucleons. 

3.2 Nonlinear equation for realistic nuclei 

Above we have considered correlated interactions of n-dipoles with a realistic (relatively dilute) 
nucleus. We have found the correlation factor Ka = (1 + ka)/ ' which is to be substituted into 
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Eq. (J2.33)) . As a result, Eq. (|2.33|) gives a generalization of the BK equation to realistic nuclei with 
correlations. As was already mentioned in the introduction, Eq. (j2.21|) respects both small and 
large limits of Ka- It is most important to stress that the correlation factor Ka is always ^> 1. 
This means the nucleus correlations increase the shadowing. In other words, the BK equation 
which assumes independent interactions underestimates shadowing corrections for realistic nuclei. 

Let us now present some numerical estimates related to Eq. (j2.21j) . The effective nucleon radius 
found in Ref. ^Hj is R% = 3.1 GeV 2 . At zero impact parameter b t = 0, is maximal. For 

the heaviest nuclei, say gold (A = 197), k,a u (0) — 0-85 while for neon (A = 20) K Ne (0) ~ 0.35. 

In the deep saturation region a solution of Eq. (J2.21|) will saturate. The saturation bound will 
not be unity but Na = 1/Ka < 1. In order to confirm our expectations Eq. (|2.21|) is solved 
numerically for gold nucleus at bt = (Fig. 0]). The initial conditions are set at x = xq = 10~ 2 . 
Consistently with our derivation we take Na(xq, r) = ka(0) Nn(xo, r) with being a dipole 
nucleon amplitude found in Ref. [TU]. The function Na shown in Fig. |U indeed saturates at the 
value 1 / Ka in accord with the expectations. 

3.3 Non- linear equation and enhanced diagrams 

The BK nonlinear evolution equation sums the 'fan' diagrams of Fig. 0-1. It is important to 
understand the limits of its applicability and develop methods to include necessary corrections. 
Two types of diagrams which are not included in the BK equation are show in Fig. ISJ-2 and 
Fig. 0-3. These are semi-enhanced diagrams including multiple Pomeron vertices (Fig. [5J-2) and 
the enhanced diagram with Pomeron loops. 

The diagram of Fig. 0-2 can be neglected in pre-saturation region only, when Na 1. Thus, 
in spite of its importance to high energy phenomenology, the validity of the BK equation is limited 
to r < 1/Qs- In the next section we will propose a new non-linear equation, which is supposed 
to sum the 'fan' diagrams as well as those of Fig. [SJ-2. 

The assumption of large ka ^ 1 was essential in arguments given in Ref. JH] i n order to 
neglect the enhanced diagrams of Fig. 0-3. Indeed, the 'fan' diagrams dominate in this limit 
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[2*4] . As we have argued above, for realistic nuclei < 1 and the argument is not applicable. 
Nevertheless we are still convinced that the Pomeron loop diagrams are subleading provided the 
projectile dipole has a small size. 

Following Refs. IH| we can estimate a contribution of the enhanced diagram of Fig. [3J-3. 
The ladders are associated with the BFKL Pomeron having the property discussed earlier: 

N oc e uoy (3.37) 

with co> oc a s . The integration over y\ and 1/2 leads to 

^enhanced 1 2o, (Y-0) ^ 2 „o F (333) 

The typical values of F — ?/i ~ l/u(r) ^> 1 while y 2 — ~ l/u(R) ~ 1. Therefore, the upper 
ladder which is in the region where the QCD coupling is small, describes a high energy process for 
which the BFKL equation can be trusted. On the other hand, the lowest ladder corresponds to the 
interaction at rather low energies where high energy asymptotic expressions are not applicable. 
Consequently, the enhanced diagram of Fig. |3J-3 should be viewed as degenerate with the 'fan' 
diagram of Fig. 0-1. A more accurate estimate of Fig. |3}3 involves integrations over dipole sizes in 
the vertices but they do not alter the conclusion (see Refs. [DI2] for details). The argument above 
depends crucially on the running QCD coupling and the fact that a s (r) <C oz s (R). Concluding we 
expect the equation summing 'fan' diagrams only be justified for nucleon targets if the projectile 
dipole r is small, such that a s (r) <C 1. 

It is worth to emphasize that even the 'fan' diagram contribution depends on the model for 
the 'two ladders-nucleus' vertex (see Fig. 0-1). In this paper we have assumed a model which 
graphically is represented in Fig. [SJ-4. This model includes the eikonal-type of interaction with the 
nucleon. 



4 Glauber rescattering and QCD generating functional 

As was discussed above, Eq. (I2.28|) for generating functional, which is equivalent to the one 
obtained by Mueller, describes the 'fan' diagrams only. The fact that this equation does not 
describe the semi-enhanced diagrams of the type shown in Fig. |3J-2 is an apparent shortcoming. 
At dipole sizes r > R s (x) = 2/Q s (x) these diagrams are of the same order as the 'fan' diagrams 
and they have to be taken into account in order to obtain a reliable description of the saturation 
domain. In this section we propose a generalization of the QCD functional approach to such type 
of interactions. To simplify the presentation we will not consider target correlations in this section. 

In order to solve the problem we have to find an effective vertex for one dipole splitting to 
many. This vertex can be deduced from the Glauber-Mueller formula fI5\ I2HJ 12Z1 which describes 
propagation of a dipole through large nucleus. It is important to stress that this formula was 
proven in QCD. In our approach, the effective vertex will be obtained as a chain of multiple 
rescatterings involving a dipole splitting to two dipoles. 
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Fig. EJ3 



Fig. EJ4 



Figure 5: The 'fan' diagrams (Fig.^jT and Fig. 0-2) for dipole-nucleus interaction. First enhanced 
diagram (Fig. 0-3). Fig. 0-4 shows a model for the interaction of the BFKL Pomeron with nucleus 
target. 
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The Glauber-Mueller formula reads 



N GM {y,r,b t ) = (l - e - NBFKL ^^) (4.39) 

with N BFKL being a solution to the linear BFKL equation. We can use Eq. ()4.39|) to define 
initial conditions for a new functional. ]\f BFKL describes a single Pomeron exchange which in the 
language of section 2 should be replace by 1 — u. Basing on the relation Eq. (j2.14|) we propose 
the following initial conditions at y = Y: 

Z{Y-y = 0, r, b t ; [u]) = e^" 1 . (4.40) 

At u — 1 , Z has the correct normalization: Z — 1. 

Eq. (14.40)) can be given the following interpretation. Contrary to the case considered in section 
2 where we had a single dipole at the beginning of the evolution, Eq. ()4.40|) describes initial state 
in which many dipoles can be created. Z represents a probability distribution for a dipole number 
in the initial state. It can be rewritten as 

ZiY = y,r,b t ; [u]) = E e ~ <n> ~ < T >} " ( 4 - 41 ) 



n=0 n - 



The average number of dipoles < n > = j-Z\ u= x = 1. Eq. (j4.41|) is a typical Poisson distribution 
with e~ <n> < n > n jn\ = e _1 jn\ being a probability to find n dipoles (of the same size) in the 
initial state. 

As in section 2, the overall probability conservation implies that at u = 1 

Z(Y-y,r,b t] [u = l]) = 1 (4.42) 

to be fulfilled at any rapidity. 

We are going now to find a new effective vertex V for one dipole splitting to many. The strategy 
is following. We consider a two step process. First, a dipole of size f decays into two dipoles of sizes 
f and f — r 1 '. This process is described by the vertex V\^i introduced in Eq. ()2.27|) . Then each 
produced dipole is considered not as a single dipole but rather having a Poisson distributed dipole 
multiplicity of equal size dipoles. This step is described using the Glauber - Mueller formula. 

The amplitude N for this two step propagation can be written in the form [27j 

= -a. W (r) (l - e- NBFKL ^) + (4.43) 



dy 

a. r r 2 d 2 r' 



2tt 



r r*dr / 



-{N BFKL (y,r') + Af SFKL (j/,r-r')} 



The first term in Eq. ()4.43|) is a possibility for the parent dipole (r) not to decay while the second 
term describes the decay. The probability for this decay is given by the square of the wave function 
1^1 2 = a s r , 2 ^_^2 ■ Two fast dipoles thus created then percolate through the target preserving 
their sizes. This process is described by Eq. 1)4.43)1 . 
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It is important to explain the difference between Eq. (|4.39|) and Eq. (|4.43|) . Eq. (|4.39|) describes 
a propagation of a dipole with the size r through target without decay. Alternative it can be 
viewed as including decays but with decay products being identical dipoles with the same size and 
rapidity On the other hand, Eq. ()4.43|) describes the two stage process. The physics of Eq. ()4.43|) 
could be clarified if we rewrite it in the following form 

dN ^ r) = -a s u(r)N GM (y,r) + (4.44) 
dy 

a, r r 2 d 2 r' 



(N GM {y, r') + N GM {y, f - f) - N GM {y, r') N GM {y, f - 7*)) 



2 7r J r' 2 (r — r 1 ) 2 

We would like to emphasize that the splitting kernel is the same LO BFKL kernel as appears in 
the BK equation. 

We are now in the position to extract from Eq. ()4.43|) an effective vertex V. We interpret the 
Glauber- Mueller formula as an annihilation of a scattering dipole and a creation of a number of 
identical new dipoles. Basing on our previous experience in section 2 we substitute in the r.h.s. of 
Eq. (J4.43|) 1 — u instead of N BFKL (see Eq. (J2.13)) and Ref.pH]) and obtain the generalized vertex 

V(r; [«]) = a s (- u(r) + / ^f/ e^' 1 + A *_ (4.45) 

\ J 2tc r' z (r — r) z J ou{r) 

The vertex V being expanded in u — 1 reproduces the vertices Vi_>i and Vi_»2 as a correct limit. 

Having obtained the vertex Eq. (|4.45j) . the linear equation for the generating functional (see 
Eq. (|2.28|l ) can be generalized straightforwardly: 



dZ 
dy 



J d 2 r' V(r'; [v]) Z. (4.46) 



As was mentioned above, if the vertex V is expanded in u — 1, Eq. (J4.46J1 reduces to Eq. (|2.28j) 
and sums the diagrams of Fig. Q 

Following the same steps as in section 2, Eq. ()4.46|) can be reduced to a non-linear equation. 
At Y — y the initial condition of Eq. (J4.4U)) lead to the relation 

S -A = Z. (4.47) 
o u 

Since Eq. (I4.4fi|) can in fact be rewritten as ordinary differential equation, we can use the relation 
(|4.47|) at any rapidity. This results in the new non-linear equation for the functional 

d In Z (Y — v, f, b t ; \u\) , . r , s 
K y, , t,l \) = _ w(r) z{Y -y, f, b t ; [u]) (4.48) 

ot s o y 

+ / ^ a( f Z(Y-y,7 7 b t ; [u]) x Z (Y - y, r - b t - [«]) . 
J 2 7i r z (r — ry 
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Eq. (J4.48j) is a new equation summing up not only the 'fan' diagrams of Fig. 0-1 but also of the 
type displayed in Fig. |SJ-2. Eq. ()4.48|) compared to Eq. ()2.29|) has an extra power of Z which we 
are not able to explain using a physical intuition. Nevertheless we believe Eq. ()4.48|) is correct in 
describing all types of 'fan' diagrams including multiple Pomeron vertices. This equation has a 
very good chance to describe the saturation region. As we have argued that inside of the saturation 
region all diagrams of semi-enhanced type are essential. 

Let us finally write down a new non-linear equation for the scattering amplitude defined in 
Eq. ()2.30|) . Using the relation Eq. ()2.14|) between the generating functional Z and the amplitude, 
the equation for the latter reads 



dN A (Y, r- b t 
a q dY 



[1 - N A (Y, r; b t )} x (4.49) 



d 2 r' r 2 



[2 N A (Y, r'; b t ) - N A (Y, r; b t ) - N A (Y, r'; b t ) x N A (Y, ?-7;b t )). 



2 7r r' 2 (r* — fy 

The initial conditions for Eq. ()4.49|) are the Glauber- Mueller formula N A (Y = y ) = N GM (y ). 
Eq. ()4.49|) has the correct linear limit (LO BFKL) when N A is small. There is a region, which 
is not defined by any expansion parameter, where N A <C 1 but N\ should not be neglected 
compared to N A . In this region Eq. (J4.49|) reduces to the BK equation. 



Quantitative estimates 

In order to estimate the impact of the newly accounted diagrams of Fig. |3J-2 we again consider 
the simple toy model of section 2.1. In this model Eq. ()4.48|) reads 

- = Z 2 (Y)(1 - Z(Y)). (4.50) 

Let define Z = Z(Y — y — 0). Eq. (|4.5()|) can be solved in the implicit form: 

l_^l e yz = Z0_^J- e !/Z e . Y (4 51) 



It can be seen that Z —>■ and iV^ = 1 — Z ^ lasF^ oo. In Fig. |Ulwe compare the amplitude 
N A calculated from the functional of Eq. ()4.51|) with the amplitude obtained from Eq. ()2.10j) the 
latter being a solution to the BK equation of this simplified model. The Fig. H shows a dramatic 
weakening of the shadowing due to the contribution of the Fig. 0-2 type diagrams. This effect 
starts to be noticeable when the amplitude N A ~ 0.5. This point is usually associated with the 
transition to the saturation region. As rapidity increases the effect strengthens and as a result, 
the amplitude approaches the black disc limit = 1 much slower compared to the standard 
behavior driven by the 'fan' diagram of Fig. 0-1. 
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Figure 6: The amplitude ob- 
tained in the toy model. The 
dotted lines show the solutions of 
the BK equation while the solid 
curves correspond to the solution 
of Eq. (j4.51|) . Four different val- 
ues for initial conditions are plot- 
ted. 



5 Results and discussions 

In this paper we developed a new approach to high density QCD (Color Glass Condensate) based 
on a linear evolution equation for QCD generating functional. We obtained three main results. 

First, we derived Eq. ()2.28j) which is a linear evolution equation for the generating functional 
summing all possible 'fan' diagrams (see Fig. EJ-l). This equation involves functional derivatives 
with respect to initial conditions which is the price paid for having a linear equation for non-linear 
dynamics. The linear equation for generating functional can be rewritten in a non-linear form 
reproducing the very same non-linear equation as derived by Mueller in Ref. |2Uj . 

The second result of this paper is a generalization of Eq. (J2.28|) . Eq. (|4.48|) is a new linear 
equation which incorporates the Glauber-Mueller rescatterings and this way sums all possible 
diagrams of the type shown in Fig. I5J-3. Eq. ()4.48|) can be reformulated as a new non-linear 
equation. The latter has a good chance to be the correct non-linear equation describing the 
parton collective phenomena such as Color Glass Condensate. 

From the equations for the generating functional we were able to obtain corresponding linear 
equations for the scattering amplitude. These equations involve functional derivatives with respect 
to initial conditions. In principal, this should not be an obstacle since the initial conditions are 
usually known and functional derivatives can be replaced by ordinary derivatives. The striking 
advantage of the linear formulation that it allows us to address a question of target correlations. 
Namely those correlations which occur in the interactions of 'wee' partons with the target. 

It follows from our analysis that the linear equation can be reformulated as standard non- 
linear evolution in two cases only. The first is when the interactions of 'wee' partons are fully 
uncorrelated. When only the 'fan' diagrams are summed this results in the BK equation. The 
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second case includes target correlations but in a most simplistic way in which only one correlation 
parameter is introduced independently of number of dipoles participating in the interaction. We 
firmly believe that within the approach based on linear equation for the scattering amplitude we 
can address a systematic analysis of target correlations and their influence on the value and energy 
dependence of the scattering amplitude. 

The third result of this paper is in the application of the ideas developed above to interactions 
with realistic nuclei. We argued that it is not correct to treat realistic (dilute) nuclei relying on 
the large A approximation. In our work we relaxed this approximation and instead introduced 
nucleus correlations in a model dependent way. The proposed model is such simple that the linear 
equation for the amplitude can be brought to a non-linear form. We obtained this way a new 
non- linear equation (Eq. 1)2.21)1 ) which is a generalization of the BK equation for realistic nuclei. 
We hope that this new equation would lead to more reliable calculations for ion-ion collisions and 
would influence our understanding of the RHIC data. 
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